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Abstract 

We report the first observation of the decay B° — ► D^D*^ with the Belle detector at the 
KEKB e + e _ collider operated at the T(4S*) resonance. The sum of branching fractions B(B° — > 
D+D*-)+B(B° -► D-D*+) is measured to be (1.17±0.26±g;||) x using the full reconstruction 
method where both charmed mesons from i?° decays are reconstructed. A consistent value ((1.48 ± 
0.38ta3i) x 10- 3 ) is obtained using a partial reconstruction technique that only uses the slow pion 
from the D*~ — > D°ir~ decay and a fully reconstructed D + to reconstruct the B°. 

PACS numbers: 13.65.+i, 13.25. Gv, 14.40. Gx 
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Mixing-induced CP- violating asymmetries have been recently observed by the B factory 
experiments using modes with charmonia in the final state such as B° — > J/ipK s Jl[|. With 
little theoretical ambiguity, the observed asymmetries determine sin20i, where 0i is one 
of the angles of the unitary triangle of the Cabibbo-Kobayashi-Maskawa matrix |2[ in the 
framework of the Standard Model. Further understanding of CP- violation and probes of 
new physics require measurements of CP- violation in additional decay modes. One class 
of such modes, the doubly-charmed decays B° — > D^D^*\ has attracted a great deal of 
attention 0]. While measurements of asymmetries in these modes are sensitive to sin20i, a 
deviation from the expected value is possible due to penguin contributions with additional 
phases from new physics. Among the doubly-charmed decays, B° — > D^D*^ is the most 
promising because neither an isospin analysis nor an angular analysis is necessary to measure 
<pi Q. So far only an upper limit for the branching fraction of B° — > D^D* T of 6.3 x 10~ 4 
was obtained by the CLEO collaboration ||, although a naive expectation for this branching 
fraction is 1.0 x 10 -3 , scaling from the well measured branching fractions for the Cabibbo- 
allowed processes B° -> D+D*~ and B° -> D~D*+ [|. 

The data used for this analysis were taken with the Belle detector [/]] at the KEKB 
asymmetric e + e~ (3.5 on 8.0 GeV) collider [§]. The integrated luminosity used for this 
analysis is 29.4 fb _1 collected at the T(4i?) resonance (-y/s = 10.58 GeV) (referred to as 
on-resonance data) and 3.0 fb _1 at a center-of-mass (CM) energy just below the threshold 
of BB pair production (referred to as continuum data). The on-resonance data correspond 
to 31.9 million BB pairs. The Belle detector is a large-solid-angle magnetic spectrometer 
that consists of a three-layer silicon vertex detector (SVD), a 50-layer central drift chamber 
(CDC), a mosaic of aerogel threshold Cerenkov counters (ACC), time-of-flight scintillation 
counters (TOF), and an array of CsI(Tl) crystals (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron flux- return located outside of the 
coil is instrumented to detect Ki mesons and to identify muons (KLM). 

In the full reconstruction method we consider the decay B° — > D + D*~ || followed by 
D*~ — > D°7r~. Five decay modes are used for D° reconstruction: K + it~, K + 7t~7t~7t + , 
K s tt + tt~ , K + 7T~7T° and K + K~; D + mesons are reconstructed through K~tt + tt + and K s ir + 
decay modes. 

Charged tracks are required to be consistent with the hypothesis that they originated from 
the interaction point in the r-<p plane. Charged kaon candidates are required to be positively 
identified based on the combined information from the ACC, TOF and CDC dE/dx systems. 
Only charged pion candidates involved in D + reconstruction are required to be positively 
identified to suppress the feed-down from the Cabibbo- allowed B° — > D^ + D*~ decays. 

Kg candidates are reconstructed from 7t + tt~ pairs with a common vertex that is displaced 
from the interaction point in the r — plane and requiring an invariant mass within ±3cr 
(±7.2MeV/c 2 ) of the nominal K s mass. 7r° meson candidates are reconstructed from pairs 
of photons, each with energy greater than 0.1 GeV, that have an invariant mass within ±3cr 
(±13.5MeV/c 2 ) of the nominal 7T° mass. The D° and D + candidates are required to have 
an invariant mass within ±2cr of the nominal masses (the mass window depends on the 
mode and varies from 8 to 22 MeV/c 2 ), where a is the channel-dependent mass resolution. 
D*~ candidates are formed from a D° and a slow pion candidate (ftJi ow )- We require the 
mass difference between D*~ and D° to be within ±llcr(~ 5MeV/c 2 ) of the known mass 
difference in order to accommodate large non-Gaussian tails in the D*~ — > D° irJ low mass 
distribution, which contain 16% of the signal. 

B° — > D + D*~ candidates are selected using the CM energy difference, AE = Eb — 



4 




FIG. 1: Kinematical distributions of B° — > D + D*~ candidates : a) scatter plot of AE versus M bc ; 

b) AE projection for M bc signal window, where the curve represents the fit described in the text; 

c) Mb c projection for AE signal window, where the curve represents the fit described in the text; d) 
cos# distributions for the data (points with error bars) and for the MC signal (shaded histogram). 



EbJami an d the beam-constrained mass M bc = J (E^^) 2 — P|, where E B (P B ) is the en- 
ergy(momentum) of the B candidate and E^ m is the CM beam energy. The scatter plot 
of AE versus M bC) and the projections onto AE and M bc are shown in Fig. [pa, 6 and c, 
respectively. If more than one candidate is found in an event, we use one of several criteria 
(e.g. best vertex-constrained fit of the D meson, number of CDC hits, best K s or 7r° fit), 
depending on the kind of ambiguity, to make a selection. For example, the best vertex- 
constrained fit is used to select among multiple D meson candidates. We define the signal 
region by \AE\ < 0.02 GeV and 5.2725 < M bc < 5.2900 GeV/c 2 , corresponding to ~ ±2cx 
and ~ ±2.5a respectively. Forty candidates are found in the signal box. Fig. [I]d shows the 
background subtracted helicity angle (6) distribution of the D*~ candidate. The distribu- 
tion follows the cos 2 9 distribution that is expected for a pseudoscalar to pseudoscalar vector 
decay. 

The signal resolution in M bc is dominated by the beam energy spread, while the AE 
resolution is dominated by the energy resolutions of D + and D*~ . Therefore the M bc distri- 
bution is less channel dependent and is more suitable for the extraction of the signal yield. 
We fit the M bc projection with a Gaussian and the ARGUS background function JTOj. The 

result is 29.6 ±6.6 signal events. The statistical significance, defined as \J— 2 ln(£(0) / £ max ) , 
is 7.0a, where £ max (£(0)) is the maximized likelihood with(without) the signal contribution. 
A consistent yield (28.0 ± 7.0) is obtained by fitting the AE distribution with a Gaussian 
function over a linear background. 

As possible sources which can contribute to the signal peak, we consider B — > D + nJ low D° 
decay, where D + and 7iJ [ow come from resonant states (Di(2420) or _D 2 (2460)) or a non- 
resonant state. The resonant state contribution is checked using the invariant mass of the 
D + and 7r^ ow . All the candidates have invariant mass below 2.2GeV/c 2 , hence there is no 



5 



TABLE I: The B reconstruction efficiencies (%) for the full reconstruction. 



Channel 




K + 7T IT 7T + 




K+tt-it K+K- 




10.7 
12.4 


3.81 
4.39 


3.88 
4.24 


2.60 7.78 
2.89 8.83 



indication of resonant states of D + and tt~i ow . The non-resonant contribution is checked 
using the Monte Carlo simulation (MC) and found to be significantly suppressed by the 
D*~ selection requirements. 

The reconstruction efficiencies for each of the ten subchannels are listed in Table | (the 
MC statistical errors on the efficiencies are negligible). For the branching fraction calculation 
B(T(AS) — ► B°B°) is assumed to be 0.5. The systematic error is dominated by the uncer- 
tainty in the charged track reconstruction efficiency (2% for each track). The uncertainty 
in the signal yield is estimated by varying the signal and background parameters in the fit. 
Contributions to the systematic error are summarized in Table |I|. Finally, the branching 
fraction of B° -> D ± D* T is determined to be (1.17 d= 0.26±g;^) x 10 3 , where the first error 
is statistical and the second is systematic. 

We also use the D*~ partial reconstruction technique previously used by the ARGUS and 
CLEO collaborations [II], [12|, [13|] which benefits from a higher B° reconstruction efficiency. 
B mesons produced in T(4S) decays have low CM momenta (~ 340MeV/c) resulting in 
charmed mesons from B° — > D + D*~ being almost back-to-back in the CM frame. The 
direction of the ft~i ow from the D*~ decay well approximates the direction of the mother 
particle due to the small energy release in the D*~ decay. Thus the CM angle a between 
tiJi 0W and D + can be employed as a signature of this decay. Another signature of the decay 
B° — ► D + D*~ is the D*~ polarization angle (8) which can be calculated in the partial 
reconstruction technique using kinematical constraints. We use the CLEO definition of the 



polarization angle in a partial reconstruction [13 



2P; 2i3 D .Y D -M D .P;' 

where P n (P*) is the k~i ow momentum in the CM frame (D*~ rest frame), and E* (E* D ) is 
the energy of the tt~ 1ow (D°) in the D*~ rest frame. Pd, Pd* and 7d* correspond to the 
D° momentum in the CM frame and the Lorentz factors for D*~ , calculated using energy 
conservation and the known B° CM energy. The polarization angle distribution should 
exhibit a cos 2 9 behavior due to the zero helicity state of the D*~ . 

To suppress the backgrounds, which are much larger than for the full reconstruction 
method, tighter D + meson selection criteria are applied. D + mesons are reconstructed 
using the K~7t + tt + decay mode only. Tracks from D + meson decays are required to have 
associated SVD hits to guarantee a valid vertex reconstruction and are refitted to a common 
vertex. The good quality of the vertex fit is required. We exploit the relatively large D + 
decay length and require the decay path direction to be consistent with the momentum 
direction of the D + : cos(Vd, Pd)t-4> > 0, where Vp is a vector from the interaction point to 
the reconstructed D + vertex. This requirement has an efficiency of 87% for real D + s and 
suppresses the combinatorial background by a factor of two. The D + CM momentum is 
required to lie in the interval 1.63 < P^+ < 1.97GeV/c and the vr^ OU) is required to have 
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FIG. 2: a) and b) cos a distributions for the leptonic sample, in regions A and B respectively; 
c) and d) cos a distributions after continuum subtraction in the non-leptonic sample. Shaded 
histograms show the expected background distributions for generic BB MC. The curves represent 
the fit described in the text. 



CM momentum smaller than 0.2 GeV/c; these requirements are the kinematic limits for the 
studied decays. 

Some of the large backgrounds exhibit a similar angular tendency to the signal, e.g. the 
cc continuum background. Continuum events are partially removed by requiring the ratio of 
the second and zeroth Fox- Wolfram moments [14] to be smaller than 0.4. To disentangle the 
signal and the continuum background, the whole data sample is divided into two subsamples: 
the sample with a high-momentum lepton in the event (referred to below as "leptonic"), and 
the sample without high- momentum leptons ( "non-leptonic" sample) . Lepton candidates are 
required to be positively identified as either a muon or an electron and to have CM momenta 
larger than 1.1 GeV/c. 

In the leptonic data sample, the presence of a high-momentum lepton in the event sig- 
nificantly reduces the continuum background. This background can be further suppressed 
using the angular correlation between the lepton and the D + candidate in the continuum. 
We require this angle to satisfy —0.8 < cos(Pd+,P^±) < 0.9. The continuum contribution is 
estimated to be smaller than 1% from the analysis of both continuum data and MC simula- 
tion. In addition, this requirement removes a significant part of the BB background where 
the lepton and the D + originate from the same B meson. 

The selected D + 7tJ loU] combinations are divided into two independent regions of the polar- 
ization angle. The region 0.5 < | cos#| < 1.05 (referred to below as region ^4) contains 85% 
of the signal events, while the region | cos#| < 0.5 (region B) is dominated by background. 
The distributions of cos a for regions A and B are shown in Fig. b. The points with error 
bars are the data while the generic BB MC background distributions are superimposed and 
shown as shaded histograms. The B° — > D^ + D*~ and B° — > D* + D*~ decay channels are 
excluded from the generic MC and are studied separately. 

In the non-leptonic sample, the continuum data distributions of cos a are scaled by a 
factor to account for the difference in the relative luminosities and cross-sections between 
on-resonance and continuum data and are then subtracted from the on-resonance data. The 
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FIG. 3: Distributions of cos 9 : a) leptonic sample; b) non-leptonic sample. Shaded histograms 
show the signal MC. 

resulting cos a distributions for regions A and B are shown in Fig. ^c,d. The data are the 
points with error bars while the generic BB MC background distributions are superimposed 
and shown as shaded histograms. 

The contamination from B° — ► D^ + D*~ and B° — > D* + D*~, which produce similar 
peaks in the cos a distribution are studied using the MC and normalized to the branching 
fractions measured in |5], [Tj|. We assume the D*~ is not polarized in these channels; the 
effect of possible D*~ polarization is included in the systematic error. The contributions of 
B° _> £)(*)+£)*- in region A are estimated to be 4.0 ±0.7 and 27 ± 5 events for leptonic 
and non-leptonic samples respectively; the B° — > D* + D*~ contributions are 3.0 ± 1.2 and 
20 ± 8 events. Several other correlated backgrounds are studied in detail. The Cabibbo- 
allowed decays B — > D*~ D^ + K result in D + CM momentum below our requirement and 
are completely rejected by the selection; for B — > Di(2A2ti)(Dl(2AQti))D^ decays the D + 
CM momentum is also shifted to lower values and these modes are strongly suppressed. The 
decays B -> L> w ~ p and B -> D^*> ~a\ are studied with different polarization hypotheses 
and are found not to produce a peak at cosa= —1. The possible contribution of these 
background sources is included in the systematic error. 

Finally, each pair of cos a distributions, Fig. @(a,&) and (c,d), is fitted simultaneously to 
signal and background functions. The signal is parameterized by an exponential function 
which is a good description of the MC signal. The ratio of signal events in regions A and B 
is fixed according to the MC. The B° — > D^ + D*~ and B° —> D* + D*~ background contri- 
butions are fixed in each region separately. The shape of other combinatorial backgrounds 
is parameterized by a second order polynomial that is assumed to be the same for regions 
A and B, in agreement with the MC. The fit finds 37.6 ± 13.2 and 283 ± 103 signal events 
in the leptonic and non-leptonic samples, respectively. As a cross-check, we fit the cos a 
distribution in eight bins of cos 9. The resulting polarization angle distributions are shown 
in Fig. |3|a and b for the leptonic and non-leptonic samples, respectively; the expected cos 2 9 
behavior can be seen in both samples. 

One of the largest contributions to the systematic error in the signal yield is the uncer- 
tainty of the background shape. We use different background parameterizations to check 
the stability of the result against the choice of the background function. We also perform a 
fit with the background shape fixed from the generic BB MC. A summary of the systematic 
error sources is given in Table ||. 

The overall efficiencies estimated for the leptonic and non-leptonic samples are equal 
to 8.4 x 10~ 4 and 5.7 x 10 -3 , respectively. We calculate the sum of branching fractions 
B(B° — > D + D*~)+B(B° — > D~D* + ) from the measurements with leptonic and non-leptonic 
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TABLE II: Systematic errors on the branching fraction (%). 



Error source 


Full 


Partial 


Reconstruction efficiency 


±17 


±12 


Particle identification 


±6 


±4 


Signal and background shape 


+1 

-7 


+13 
-15 


Correlated background 




+2 
-3 


Number of BB pairs 


±1 


±1 


D branching fractions 


±7 


±7 


Total 


+19 
-21 


+19 
-21 



samples to be (1.41 ± 0.52±g;§2) x !0~ 3 and (1.57±0.57lg;^) x 10 -3 , respectively, where the 
first error is statistical and the second is systematic. Averaging these two measurements we 
calculate B(B° -> D + D*-)+B(B° -> D~D* + )=(1A8 ± 0.38±g;§f) x 10" 3 . 

In summary, we report a measurement of the sum of the branching fractions for the decays 
B° — > D + D*~ and B° — > D~D* + using two analysis methods. The sum of the branching 
fractions B(B° -> D + D*-)+B(B° -> D~D* + ) is calculated to be (1.17 d= 0.26±g;l§) x 10~ 3 
using the full reconstruction method and (1.48 ± 0.38lo;li) x 10~ 3 using the partial recon- 
struction technique. For calculations of the average branching fraction in this channel, the 
result from full reconstruction should be used. The partial reconstruction result is a consis- 
tency check and the sample obtained using this method will be used to improve statistics 
in future measurements of CP-violation. Our results are higher than the upper limit from 
CLEO; they are however consistent with the expectation based on B° — > D*D^*\ 
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